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I .  Introduction 

This  final  report  summarizes  the  research  conducted  under 
the  Contract  DAAG  29-79-C-0075  during  the  period  July  1,  1979 
to  June  30 ,  1980.  It  was  the  goal  of  this  research  to  in¬ 
vestigate  short,  intense  submillimeter  pulse  generation  ob¬ 
tained  by  mode  locked  optical  pumping.  Efforts  were  made 
towards  demonstrating  short  pulse  generation  at  1.2mni  and 
towards  the  general  characterization  of  such  an  approach  in 
experiment  and  in  theory;  the  latter  two  of  which  will  be 
outlined  in  the  following  two  sections.  The  last  section  is 
devoted  to  the  summary  of  this  report. 

II .  Experiment 

The  technique  that  has  been  explored  in  this  research 

was  in  the  use  of  mode-locked  optical  pumping  of  low  pressure 

gases,  which  has  been  previously  demonstrated  successfully  for 

subnanosecond  pulse  generation  at  385  yin  in  O2O  and  496  pm  in 
12  1 

C  H^F.  Recently,  this  approach  has  been  extended  to  other 

transitions,  e.g.,  66  ym  in  D2O,  193  ym  in  C^^H^F  and  151  ym 
2 

in  NH^.  This  study  was  aimed  at  not  only  the  generation  of 
short  intense  pulses  at  1.2mm  in  C^H^F  but  also  the  mathema¬ 
tical  characterization  of  such  an  optical  pumping  approach. 

The  experiment  was  performed  in  a  simple  single-pass 

configuration  as  is  shown  in  Figure  1.  Here  the  hybrid  TEA- 

3  4 

CO 2  laser,  which  was  used  for  other  experiments,  ' 


was 


FIR  CELL 
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converted  to  a  mode-locked  laser  with  the  insertion  of  a  Ge 
acousto-optic  modulator.  This  laser  provided  a  very  stable 
train  of  mode-locked  pumping  pulses,  whose  pulsewidth  was 
variable  from  2  nsec  to  5  nsec,  depending  on  the  RF  power 
applied  to  the  modulator. ^  The  Pyrex  FIR  waveguide  was  term¬ 
inated  with  Si-Brewster  window-couplers  at  both  ends.  The 
path  of  the  C02  laser  beam  was  such  that  it  entered  the  FIR 
cell  via  the  input  coupler,  propagated  to  the  end  and  finally 
exited  through  the  output  coupler.  Of  the  two  possible  propa¬ 
gation  components  of  the  FIR  wave,  the  backward  wave  was  fed 
back  to  the  cell  by  a  reflector  at  the  entrance  end  to  enhance 
the  signal  when  the  FIR  output  was  not  appreciable,  but  other¬ 
wise  it  was  absorbed  by  an  absorber.  The  C02  laser  exiting  the 
cell  was  focused  onto  a  Ge  photon  drag  detector  to  monitor  the 
pumping  pulses.  The  FIR  pulses  were  detected  by  a  tungsten 
wire-on-nickel  MOM  diode,  and  observed  with  a  Tektronix  485 

oscilloscope  (350  MHz  bandwidth) . 

13 

The  C  isotope  of  methyl  fluoride  has  a  few  lasing  tran¬ 
sitions  at  millimeter  wavelengths  (1.222mm,  1.207mm,  1.006mm 
and  0.862mm) ®  and  at  submillimeter  wavelengths  (388  ym,  412 
ym)  due  to  the  hot-band  transition.7  Figure  2  shows  the  par¬ 
tial  energy  level  diagram  and  some  of  the  relevant  lasing 
transitions,  however  since  the  operating  pressure-  was  around 
one  torr,  the  dominant  transition  was  expected  to  be  mainly 
at  1 • 2mm. ^ 


The  pulse  width  observed  from  the  mode-locked  pumping 
scheme  ranged  from  5  nsec  to  8  nsec,  some  typical  pulses  of 
which  are  shown  in  Figure  3.  The  pulse  shape  was  determined 
to  be  very  close  to  a  hyperbolic  secant  waveform.  The  FIR 
pulse  envelope  delay  with  respect  to  that  of  pumping  pulses 
was  also  measured  at  different  pressures,  and  Figure  4  shows 
the  combined  signals  due  to  the  pump  and  the  FIR  pulses. 

This  result  is  compared  with  the  data  obtained  from  the 
swept-gain  superradiance  experiment,  which  is  shown  in  Fig¬ 
ure  5.  It  is  noticed  that  there  is  a  good  agreement  between 
the  two  experiments  suggesting  a  superradiant  like  behavior 
in  one  case.  The  pulsewidth  variation  with  the  cell  pressure 
was  not  significant  enough  to  be  resolvable  by  the  present 
detection  system.  The  scaling  law  of  pulsewidth  in  terms  of 
such  parameters  as  cell  pressure,  pump  intensity  or  energy 
and  cell  length  is  yet  to  be  confirmed  experimentally,  but 
the  following  theoretical  consideration  suggests  that  some 
pulse  width  reduction  be  possible  through  various  means.  In 
summary,  pulses  at  about  the  kW  level  and  of  duration  ^  Tj 
have  been  obtained  at  1.22mm  by  this  technique. 

III.  Theory 

In  addition  to  the  experimental  investigation,  a  theore¬ 
tical  study  was  also  considered  in  order  to  understand  the 
conversion  dynamics  of  the  mode-locked  pumping  of  low  pressure 


This  was  taken  at  1  torr 
long.  Horizontal  scale 
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Figure  4 
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FIR  envelope  delays  at  three  different  pressures; 
a)  0.1  torr,  b)  0.07  torr,  c)  0.05  torr.  The 
pressure  was  measured  with  a  capacitive  manometer 
Horizontal  scale,  0.1  ysec/div. 
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gases  and  to  assess  the  practicability  of  such  an  approach. 

In  view  of  the  fact  that  the  pulses  are  shorter  than  the  'Tj' 

time  scale  of  the  medium,  the  conversion  process  belongs  in 

the  realm  of  a  true  transient  phenomena.  There  are  a  few 

other  studies  that  have  been  considered  and  are  related  to  the 

present  situation.  One  is  the  transient  stimulated  Raman  scat- 
8-io 

tering  (TSRS)  and  the  other  is  the  swept-gain  superradiance. 

In  this  study,  the  near- resonance  nature  of  the  interaction 
precludes  the  former  from  being  directly  applied.  The  latter 
may  be  relevant  at  least  for  non-overlapping  pump  and  FIR,  but 
the  FIR  pulse  evolution  during  the  pumping  period  has  to  be  in¬ 
cluded.  This  has  been  studied  recently  in  terms  of  the  initial 

11  12 

tipping  angle  evolution  during  the  pump  pulse  period.  ' 

In  this  study  we  modelled  the  situation  as  a  resonant  in¬ 
teraction  of  two  waves  with  a  three-level  system,  and  the 
transient  analysis  was  comprised  of  solutions  of  the  semiclas- 
sical  model  based  on  the  equations  of  motion  for  the  density 
matrix  elements  along  with  the  wave  equation  for  the  FIR 
field.  The  energy  level  configuration  for  which  the  analysis 
was  based  is  shown  in  Figure  6.  After  some  standard  algebraic 
manipulations,  the  following  normalized  set  of  coupled  equa¬ 
tions  is  obtained  where  t  is  the  retarded  time,  t  »  t  -  •£-  . 

s 

Diagonal  Element  or  Population; 


i 

at  *  1  epr13 


(Nj-nJ) 
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3Es 

az  "g  R32  "  kEs 

In  these  equations,  the  following  definitions  are  intro¬ 
duced: 

Rabi  frequency:  Ep  =  uIR*£p/ft,  Eg  = 

PIR  *  IR  transition  dipole  moment  0.1D),  PpER  =  FIR  transi¬ 
tion  moment  ('v.  ID) ,  N®  *  equilibrium  population  in  the  ground 
2 

state,  g  =  and  k  ■  linear  field  loss  coefficient 

representing,  for  example,  diffraction  losses.  In  this  simpli¬ 
fied  treatment,  pump  depletion  is  neglected. 


12 


We  were  interested  in  the  solutions  of  these  equations 

in  two  different  limits.  One  was  the  time  response  of  the 

system  to  the  repetitive  pumping  by  the  periodic  pulses,  and 

the  other  was  the  FIR  pulse  evolution  both  in  time  and  space 

in  a  single  pulse  scheme.  First,  a  numerical  integration  was 

performed  on  these  equations  for  a  train  of  mode-locked  pumping 

pulses  with  a  parabolic  power  spectrum.  Here,  the  M-degeneracy 

factor  was  included  in  the  analysis  and  the  FIR  field  was 

assumed  to  have  reached  the  so-called  'steady-state',  i.e., 

space  independent.  This  limit  was  described  in  more  detail  in 
13 

an  early  report  so  only  some  of  the  salient  features  will  be 
discussed  here. 

Figure  7  depicts  typical  pump  and  FIR  pulses  (from  the 
numerical  selection)  in  the  initial  portion  of  the  train  sharing 
the  process  of  FIR  pulse  stabilization.  Once  the  FIR  has  been 
stabilized,  the  periodic  regeneration  is  maintained.  In  Fig¬ 
ure  8,  the  effect  of  the  peak  pump  intensity  variation  on  the 
FIR  pulse  shape  is  shown.  For  a  weak  pump,  the  FIR  pulse 
shape  is  closely  given  by  a  hyperbolic  secant  waveform.  As 
the  pump  intensity  is  increased,  the  delay  between  the  two 
peaks  decreases  as  does  the  FIR  pulsewidth,  and  for  even  a 

stronger  pump  the  FIR  pulse  leads  the  pump  pulse.  The  latter 

2 

has  been  observed  in  the  other  experiments.  The  main  point 
to  note  from  these  is  that  pulse  widths  of  about  half  the  pump 
pulse  width  are  predicted. 


FIELD  (sec' 


+  + 


+  + 


IO#h+  + 


TIME  (nsec) 


Figure  7.  Normalized  fields  (Rabi  frequencies)  versus  time. 

The  calculation  was  started  at  t«-3  nsec  with  an 
initial  condition  of  E  (o)-l.  For  t  >  50  nsec,  the 
system  has  stabilized.  The  pump  field  is  approxi¬ 
mately  1  MW/ cm2,  1.2  nsec  duration. 


RELATIVE  INTENSITY 

•  • 


.8 

6 

4 


49  50  51  52  53  54 

r(TIME) 


Figure  8.  Comparison  of  the  Stokes  pulses  at  different  pump 

intensities.  The  peak  pump  intensity  for  each  case 
is,  a)  105  M/cm2,  b)  10®  W/cm2,  and  c)  107  W/cm2. 
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In  a  second  inquiry,  the  analytical  solution  of  the  equa¬ 
tions  was  pursued.  Here  the  Stark  shift  term  due  to  the  FIR 
field  was  neglected  and  the  solution  was  obtained  for  the  FIR 
pulse  where  the  pump  pulse  width  was  restricted  to  t  <<  Tj. 

The  calculational  details  are  presented  in  the  Appendix.  After 
solving  the  equations  in  a  self-consistent  manner,  the  FIR 
field  is  obtained  as  Equation  (A-10) 


E_  (t,  z)  *  E  (t, o)  +  Jqz  sin 

S  S  e 

t  E  (t-t1  ,o)  I.  f  2/gzt  (t ' )  I 

— ^  si 


0  Sin  2 


/T(t') 


.  2  e(t') 

sin  — ^ —  dt 1 


where 

e  (t) 


t 

E  (t')dt',  the  pump  pulse  area, 
0  p 


Es ( t, o)  =  FIR  initiating  field  at  z=0  s  EsQ 

1^  =  modified  Bessel  function  of  the  first  kind,  and 

x(t')  *  f  sin2  9 ■■  dt". 

'0 


Invoking  the  asymptotic  limit  of  1^  for  a  small  argument, 
the  FIR  field  is  simplified  as 


EJt,z)  s  E_  +  gzE^  sin 


ill l 


a 


so 


so 


in  2mn 


sin 


sin 


Tt=tT  dt’ 


*A 
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which  shows  a  linear  growth  in  z.  For  large  agrument, 
I^(x)  -*■  (2irx)  ^e*  so  that 


sin 


e^ti 


sin 

3HEEI 


2/gzx (t ’ ) 


sin- 


[T(t)] 


77T 


dt1 


Since  x(t)  is  an  integral  of  a  squared  function  and  is  expected 
to  be  a  smooth  function  of  time,  this  result  implies  a  FIR 
field  growth  scaling  as 


Eg ( t,  z)  <*  Ego  exp 


2/gzr 


From  these  asymtotic  solutions,  we  can  see  that  the  FIR  field 

growth  is  linear  in  z  in  the  beginning  but  as  it  propagates, 

its  growth  evolves  to  exp^2/gzx j .  For  a  rectangular  pump  pulse 

with  small  area,  the  FIR  pulse  can  be  shown  to  grow  as 

exp(i/Gzt  )  where  GT2  *  9E*T2  and  is  interpreted  as  the 

steady-state  Raman  gain.  As  was  pointed  out  in  the  early 

report,  this  characteristic  implies  the  temporal  contraction 

of  Stokes  pulse  as  it  propagates.  The  sine  factor  in  the 

solution  results  in  structure  as  the  FIR  pulse,  which  is  the 

manifestation  of  the  nutation  effect  caused  by  the  pump.  This 

is  analogous  to  the  effect  of  the  material  saturation  in  the 
10 


TSRS. 


IV .  Summary 


In  concluding,  this  work  extends  the  mode- locked  optical 

pumping  technique  into  the  mm  range  with  the  production  of  ^ 

13 

5  nsec  pulses  from  C  H^F.  The  theoretical  description  of 
the  conversion  process  was  also  investigated  in  terms  of  the 
small  signal  and  large  signal  growth  of  the  FIR  field  with 
the  main  result  being  the  conversion  is  more  closely  related 
to  swept  gain  superradiance  than  to  stimulated  Raman  emission. 
With  these,  we  can  set  rough  bounds  on  the  shortest  possible 
FIR  pulse  widths  Tp  as  1.8  T2<p/Gp  <  xp  ^  tp/2. 
limit  comes  from  assuming  a  delta  function  pump  and  is  the 
steady  state  superradiance  result  while  the  upper  limit  is 
from  the  numerical  solutions  discussed  earlier.  Using  values 
appropriate  to  the  experiments,  FIR  pulses  in  the  range  of 
200-400  psec  may  be  anticipated.  What  remains  to  be  deter¬ 
mined  is  the  spatial  conversion  distance,  overall  conversion 
efficiency  and  the  pulse  to  pulse  (in  the  train)  stability  of 
the  FIR.  The  latter  is  of  interest  for  possible  encoding 
schemes . 
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V.  Appendix 

After  neglecting  the  Stark  effect  terms  due  to  the  FIR 
field,  the  density  matrix  equations  can  be  reduced  to  the  fol 
lowing  set  of  equations  appropriate  to  times  t  <<  Tj : 


3W 


—  =  E  R 


3t 


pl3 


3W 


32 


at 


-  i  Via 


3R 


13  _ 


3t 


3R 


32 


3t 


-  EsW13 


EsW32  +  2  EpR12 


(A-l) 


3R. 


12  _ 
3 1 


1  EpR32  +  \  EsR13 


where  ”  Ni‘"N3'  W32  *  N3-N2  are  P°Pulati°n  differences. 
Solving  these  equations  simultaneously,  the  polarization  ele¬ 
ment  R^2  for  the  FIR  field  is  found  to  be  given  by 


R32(t,z)  *  -  sin 


Eg  (t*  ,  Z) 


sin 


mu 

2 


dt' 


(A-2) 


t» 

where  8(t')  =  E  dt,  the  pump  pulse  area. 

1  o  p 

Substituting  this  into  the  wave  equation  for  the  FIR 
field  and  neglecting  the  loss  term,  one  obtains 


r 


3Es(t,z) 


g  sin  Eg(t ' , z)  sin  dt' 


for  the  field  growth. 


>Mt,z)  = 


Eg(t ' ,  z)  sin  —  dt' 


(A- 3) 


(A-4) 


Then  iHt,z)  satisfies  the  following  hyperbolic  equation. 


.  g  sin2  iMt,z> 


The  transformation  of  independent  variable 


(A- 5) 


2  m 


sin  --^2  dt 


(A-6) 


reduces  Equation  (A- 5)  into  the  standard  form  as 


32»Ht,z] 


(A-7) 


The  boundary  conditions  are  such  that 


(t,z) 

3z 


<Mt,Z) 


(■t 

Eg(t',o)sin  - ^  dt' 


(A-8) 
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Equation  (A- 7)  can  be  solved  by  the  Riemann's  method; 
the  Riemann  function  of  which  is  found  to  be 


RU,n;t,z)  =  i0|2/g(x-5)  (z-n)] 


where  IQ  is  the  modified  Bessel  function  of  zeroth  order. 

Therefore,  the  solution  to  Equation  (A- 7)  can  be  obtained 


as 


iMt,z)  =  y(T,o)  + 


3Ift  2/gxz 

^(t-x,o)  - ^ - L  dx 


(A- 9) 


Taking  the  derivatives  of  both  sides,  from  Equation  (A-4) 
the  FIR  field  is  found  to  be  generated  by 


Eg(t,z)  =  Eg(t,o)  +  Jqz  sin 


tEg(t-t',o)  ^ 

2/gzr  ( t ' )  ) 

“TeTt^FT 
i  0  Sin  2 

/T(tri 

’  dt ' 


(A-10 ) 


which  can  be  integrated  numerically  for  non-simple  pump  pulse 
areas. 
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